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Abstract 

Although it is known that correct dynamical equations of motion for 
a nonholonomic system cannot be obtained from a Lagrangean that has 
been augmented with a sum of the nonholonomic constraint equations 
weighted with multipliers, previous publications suggest otherwise. An 
example has been proposed in support of augmentation and purportedly 
demonstrates that an accepted method fails to produce correct equations 
of motion whereas augmentation leads to correct equations; this paper 
shows that in fact the opposite is true. The correct equations, previ- 
ously discounted on the basis of a flawed application of the Newton-Euler 
method, are verified by using Kane’s method and a new approach to de- 
termining the directions of constraint forces. A correct application of the 
Newton-Euler method reproduces valid equations. 


Introduction 

Dealing with nonholonomic constraint equations within the framework of varia- 
tional methods is a controversial subject. For example, Ray [1] modifies Hamil- 
ton’s principle and augments the Lagrangean by adjoining a sum of nonholo- 
nomic constraint equations weighted with multipliers. Later, Ray [2] reverses 
himself. In the erratum he compares the correct way of dealing with constraint 
equations that are linear in the time derivatives of the generalized coordinates 
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to the incorrect approach of augmenting the Lagrangean that gives the wrong 
results, even when the constraint equations are linear. Saletan and Cromer [3] 
follow Ray and show the augmented Lagrangean gives correct equations of mo- 
tion when the constraint equations are holonomic. They conclude that no such 
augmented Lagrangean exists in the nonholonomic case, in part because they 
say that there is no way to determine initial conditions needed for the integration 
of differential equations governing the multipliers. Rosenberg [4] presents the 
same demonstration as Ray’s erratum and concludes that, although Hamilton’s 
principle may be regarded as a variational principle for conservative holonomic 
systems, it cannot be so regarded for nonholonomic systems. In an effort to elim- 
inate constraint violations, Rosen and Edelstein [5] make a proposal similar to 
that of Ray 30 years earlier; they account for nonlinear nonholonomic constraint 
equations in the same way that they do holonomic constraint equations. Hage- 
dorn [6] points out that although their approach is justified in the holonomic 
case, it is incorrect for nonholonomic constraint equations, even when they are 
linear. He demonstrates this with an example and gives the well-known result for 
the correct way to handle linear equations, which does not come from modifying 
the Lagrangean. According to Hagedorn the mistake has been repeated many 
times over the past century and the pitfall has received attention in Refs. [7], 
[8], and [9]. More recently, Flannery [10] examines the problems encountered 
by Ray and others and, after in-depth analysis, concludes “General [nonlinear] 
nonholonomic constraints are completely outside the scope of even the most 
fundamental principle of D’Alembert. The generalization of any principle based 
on [D’Alembert’s] to general nonholonomic constraints is without foundation.” 
In an earlier paper, Polyakhov [11] reaches the same conclusion: “Thus, the 
D’Alembert-Lagrange principle as well as the Lagrange ideality condition can 
be derived from Newton’s equations only for holonomic constraints. Therefore, 
the application of these results to nonholonomic systems is not justified.” 

In their response to his comments, Rosen and Edelstein offer a counterex- 
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ample purportedly showing that the approach advocated by Hagedorn leads 
to incorrect results. Their conclusion is based on a flawed application of the 
Newton-Euler method. The purpose of this paper is to demonstrate the validity 
of the approach taken by the latter author and bring to light mistakes made by 
the former authors in their development of the counterexample. 


Problems With the Counterexample 

The planar system featured in the counterexample is shown in Fig. 1. Two 
perpendicular unit vectors hi and h 2 are fixed in an inertial reference frame 
N. A smooth rod B whose axis is parallel to unit vector bi is inclined at a 
constant angle a to hi; B is permitted to translate along an axis parallel to 
112 . A particle P of mass m moves along P, and the mass of B is negligible in 
comparison to m. It is said that no forces are exerted on P other than those 
necessary to prevent it from loosing contact with B. 

Analysis is facilitated by working with two generalized coordinates q\ and 
q 2 shown in Fig. 1, where q\ is the displacement in a prismatic joint connecting 
B to N , and where q 2 is the displacement of P along the rod. Two generalized 
speeds [12] are introduced simply as u r = q r (r = 1, 2). A motion constraint is 
to be imposed upon the velocity of P in P, expressed by the relationship 

e C a u 2 ~ qi = 0 (1) 

where e is a positive constant and C a denotes cos a. The equivalence of this 
expression and the constraint equation in Ref. [6] is demonstrated presently. 
Now, the velocity N v B in N of every point B fixed in B is given by N v B = 
iqh 2 , and the velocity B w p of P in B is given by p v p = u 2 bi. Henceforth, B 
is taken to be the point of B that is coincident with P, and the velocity of P in 
N is simply N v p = N v B + p v p . The nonholonomic constraint equation (1) 
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Figure 1: A Particle Moving on a Sliding Inclined Rod 
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can thus be written in vector form as 

(iY v F_ iV v B } . 6i _^ =0 (2) 

There also exists a configuration constraint that prevents P from moving in B 
in the direction of b 2 ; at the velocity level, the holonomic constraint equation 
is expressed as 

(*v p - V).b 2 = 0 (3) 

On the basis of physical reasoning it is clear that a force normal to the rod, in 
other words a force of unknown magnitude and parallel to b 2 , must be applied 
to P in order to satisfy Eq. (3) and keep P from losing contact with the rod. In 
view of the law of action and reaction, a force of equal magnitude and opposite 
direction is thus applied to B at B. These conclusions can be reached also by a 
close examination of the form of Eq. (3), which is now undertaken so that the 
results may be used with Eq. (2) and another constraint equation to follow. 

The vector form of Eq. (3) can be used to determine directions and points 
of application of forces necessary to impose the constraint. Upon differentiation 
with respect to time in TV, one obtains 

(*a p - N a p ) • b 2 = 0 (4) 


where N a p and N a. B are, respectively, the accelerations of P and B in TV. 
(The time derivative of b 2 in TV vanishes because B has no angular velocity 
in TV.) Each acceleration is related by Newton’s second law to the resultant of 
all applied contact and distance forces. Let f p and fg denote the respective 
resultants of all contact and distance forces applied to P and B when their 
motion is not restricted by the constraint expressed in Eq. (4). In addition, 
let C p and Cg represent constraint forces that must be applied to P and B 
respectively in order to impose the constraint. By virtue of Newton’s second 
law, Eq. (4) can be rewritten as 


fp + Cp • b 2 + . (_b 2 ) = o 


m 


m B 


(5) 
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where m g is the mass of the particle at P. Any component of Cp that is 
perpendicular to b>2 will not play a part in Eq. ( 5 ); therefore, all that is necessary 
is for the constraint force Cp to be parallel to b>2- Likewise, Cp need only be 
parallel to — t>2. The first condition is expressed with the relationship 

C2 = fl 2^2 ( 6 ) 

where C2 is the constraint force applied to P and associated with the second 
of two constraint equations [Eq. ( 3 )], and where /i 2 is an as yet unknown scalar 
multiplier. Furthermore, the constraint force applied to B is /r 2 (—£>2) = — C2. 
One may simply inspect Eq. ( 3 ) to make such determinations; because N w p 
appears in a dot product with £>2, it can be said that the constraint force parallel 
to £>2, C2 = /i2t>2, must be applied to P. Likewise, because N v B appears in a 
dot product with — £>2, it can be said that the constraint force — C2 = — ^£>2 
must be applied to B at B. As has been said, these results are in line with 
physical reasoning. In Ref. [ 6 ] the magnitude of the reaction force C2 normal 
to bi is denoted by N rather than /i 2 . 

Using the same reasoning, one may inspect Eq. ( 2 ) and conclude that P 
must be subject to a constraint force Ci that is parallel to £q, 

Ci = 11 ibx ( 7 ) 

whereas a force — Ci is applied to B at B. In this case the subscript of Ci and 
fi | indicates they are associated with the first of two constraint equations. In 
practice the set of forces Ci and — Ci could be applied with a motorized gear 
attached to P moving on a track of gear teeth fixed in P; evidently the rod 
cannot be perfectly smooth as hypothesized in the problem statement, if the 
nonholonomic constraint equation ( 1 ) is to be satisfied. 

After forming the acceleration of P in A as N a p = ui 112 + t^bi, one is 
in a position to use Kane’s method [ 12 ] and form two equations of motion for 
the system S composed of P and P, F r + F* = 0 (r = 1 , 2 ). The holonomic 


December 7 , 2006 


6 


Roithmayr 


An Argument Against Augmenting the Lagrangean for Nonholonomic Systems 


generalized active forces are given by 


F r = N v r p ■ (Ct + C 2 ) + N v r B • (-Ci - C 2 ) (r = 1, 2) (8) 


and the holonomic generalized inertia forces F* are constructed according to 


so that the dynamical equations of motion for S' in A" are found to be 


where S a denotes sin a. The nonholonomic constraint force Ci contributes to 
the holonomic generalized active forces, whereas the holonomic constraint force 
C 2 does not. A third equation is needed to solve for the three unknowns tq, tq, 
and /q; it is provided by the nonholonomic constraint equation (1) expressed 
at the acceleration level, 


It is worth noting that Eqs. (7), (15), and (16) together contradict the statement 
in Ref. [6] preceding Eq. (16) therein. The acceleration of P along the path (the 
rod) is in general nonzero if the proposed nonholonomic constraint equation is 
to be satisfied; it vanishes only in the special case when the rod is stationary 

K =0). 




(r = 1, 2) 


( 9 ) 


The required holonomic partial velocities are 



ra(iq + S a ii 2 ) = 0 
m(S a ui + u 2 ) = /q 


(11) 

( 12 ) 



( 13 ) 


An analytical solution is then available, 


ill = iq 


tana 


(14) 

(15) 


e 



mC a 


(16) 


Hi = tii 


e 
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The differential equation (14) yields a closed form solution 

Ul = K x e~ pt (17) 

where p = (tana)/e as defined in Ref. [6], and the constant of integration K\ is 
the value of u\ at t = 0. This solution in turn facilitates integration of Eq. (15), 
yielding 

U 2 = ~e~ pt + K 2 (18) 

da 

where K 2 is determined once the value of U 2 at t = 0 is specified. Integration of 
the two kinematical differential equations q r — u r (r = 1,2) produces solutions 
for the generalized coordinates. 

qi = -—e- pt + K z , q 2 = ^- e ~ pt +K 2 t + K 4 (19) 

P Pd a 

where the constants of integration K 3 and K 4 can be evaluated on the basis of 
the initial conditions of q\ and q\ 2 . It must be noted that the initial values of U 2 
and qi have to satisfy Eq. (1); a similar recognition appears in Ref. [6]. 

One is now in a position to show that Eqs. (17)-(19) verify the results 
attributed in Ref. [6] to Hagedorn’s approach. First, relationships between the 
Cartesian coordinates x and y and the generalized coordinates q\ and q> 2 are 
established. 

X = CaQ2 5 x = C a U2, x = C a U 2 (20) 

y = qi+Saq 2 , y = Ui+SaU 2 , y = U! + SaU 2 (21) 

Appropriate substitution from these relationships shows that the original form 
of the nonholonomic constraint equation given in Ref. [6], y — xta,na — ex = 0, 
gives way to Eq. (1). Furthermore, Eqs. (12) and (14a) in Sec. 3 of Ref. [6] are 
recovered from Eqs. (19) here. 

X = C a q 2 = -Tl _ e -pt + K 2 c a t + K 4 c a = D 3 e~ pt +D 4 t + D 5 (22) 
p tana 

y = qi+ s a q 2 = e~ pt +K 3 + ^ e~ pt + K 2 S a t + K 4 S a =D 4 t + D 2 (23) 
P P 
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It is clear that because of the constraint described by Eq. (1), only three of the 
constants of integration Ah, K 2 , and K 4 are independent, and only three 
of the five constants Di , . . . , D$ used in Ref. [6] are independent. 

The second order differential equations (11a) and (lib) in Sec. 3 of Ref. [6] 
can also be recovered from Eqs. (14) - (16) here. Dealing with the differential 
equation for y is straightforward, 

... ^ . tana u\ . 

y = Ml + S a u 2 = Ui + S a —— = 0 (24) 

e e o Q 

The differential equation for x can be rewritten as 

X = C a u 2 = C a ^r = Ay (25) 

e C a mC a 

A relationship between the multiplier yiq used here and the multiplier A used in 
Ref. [6] is required, and can be obtained by rewriting the original nonholonomic 
constraint equation as 

— ex + y — xtan a = N v p • (— ehi) + y — xtan a = 0 (26) 

Inspection of this equation indicates that a constraint force parallel to the vector 
— ehi must be applied to P. Consequently, 

Ci = — Aeni (27) 

The projection of onto tq must be the same as that of Ci; therefore, 

Cx • bi = Ci • bi = /ii = -A eC a (28) 

Hence, Eq. (25) is rewritten 

x + — = 0 (29) 

rn 

in agreement with Eq. (11a) of Ref. [6] when m is taken as unity. 

Rosen and Edelstein reject the preceding differential equations for x and y, 
and the closed form solutions, on the basis of their results obtained with the 
Newton-Euler method. With the analysis already performed here it is evident 
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that their application of the method is flawed, and the point in their development 
where the mistake was made can be identified immediately. In what follows, a 
correct application of the Newton-Euler approach is shown to yield the foregoing 
results. 

In their Eqs. (2a) and (2b), Rosen and Edelstein do not account for the con- 
straint force needed to ensure satisfaction of their nonholonomic constraint 
equation; they only consider C2 required to bring about the configuration con- 
straint. Upon writing + C2 = m N a p , it is seen that Eqs. (2) should be 
stated 


mx = —NS a — Ae, 

my = NC a 

(30) 

= —NSa + Ml/Ccn 

m(iii + S a u 2 ) = NC a 

(31) 


Application of Newton’s second law to P must be accompanied by its application 
to B. In addition to the reaction forces — and — C2 acting at 5, a reaction 
force /13111 is applied at the prismatic joint, therefore we write fi^rvi ~ — 

C2 = rriB N a B ; however, the mass of B is neglected in comparison to m so 
/i3iii — C[ — C2 = 0 . That is, 

/i 3 + Ae + NS a = 0, —NC a = 0 (32) 


The first of these can be used if fi 3 is of interest but the relationship is not 
important in what remains to be done. The second of these reveals that N = 0 
(so fi 2 = 0) and, as an immediate consequence, y = 0, in agreement with 
what has been previously shown. Equations (31) can now be simplified, and the 
nonholonomic constraint equation (13) at the acceleration level is again brought 
to bear so that 

mC a U2 — = 0, iii + S a ii2 = 0, h 2 = 0 (33) 

constitute three equations in three unknowns, hi, and fii. They lead imme- 
diately to Eqs. (14)-(16). 
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Conclusion 

A straightforward application of Kane’s method for simple nonholonomic sys- 
tems, together with identification of the constraint forces needed to impose 
a motion constraint and a configuration constraint, are used to verify results 
obtained with what is called the regular variational approach, brought to the 
reader’s attention by Hagedorn in Ref. [6]. Further, a correct application of the 
Newton-Euler method reproduces those results. The conclusion by Hagedorn, 
Ray, Flannery, and others is thus affirmed; namely, the Lagrangean cannot be 
augmented by the sum of nonholonomic constraint equations weighted with 
multipliers, regardless of whether or not such equations are linear in the time 
derivatives of the generalized coordinates. 
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